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■ Sound attenuation and internal friction coefficients are calculated for a real- 

■ istic model of amorphous silicon. It is found that, contrary to previous views, 
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thermal vibrations can induce sound attenuation at ultrasonic and hypersonic 
frequencies that is of the same order or even larger than in crystals. The rea- 
j-^ ■ son is the internal-strain induced anomalously large Griineisen parameters of 

o ■ 

O ■ the low- frequency resonant modes. 
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Sound attenuation in glasses is poorly understood. This is because many competing 
factors lead to sound wave damping. Most important are thermally activated structural re- 
laxation, hypothetical tunneling states, topological defects, and thermal vibrations. Sorting 
out different contributions for a given temperature T and sound wave frequency z/ = VL/2n 
is a difficult task. 

Experiments show the following features: (i) At temperatures T ^ 10 K and ultrasonic 
frequencies (10 MHz to 1 GHz) the sound attenuation coefficient r(T) exhibits a small, 
frequency-dependent peak |]l|. (ii) At higher temperatures, between 10 and 200 K, another 
peak in T{T) develops whose center increases only moderately when v increases. The peak 
broadens at hypersonic frequencies and is not seen above 100 GHz [^3. As a function 
of frequency T{u) ~ z/ at the peak temperatures (iii) At hypersonic frequencies T{T) 
appears to be almost independent of (or slightly increasing with) T above the peak (ii) to at 
least 300 K ||JD. (iv) Room temperature T{v) ~ z/^ from at least 200 MHz [|[; this depen- 
dence continues for up to about 300 GHz and seems valid for any temperature above 
the peak (ii). Finally, (v) the attenuation coefficients for longitudinal (F^) and transverse 
(Ft) waves are similar. 

While the low-temperature behavior (i) of F is understood based on the interaction of 
sound waves with tunneling states [|I|, features (ii) through (v) lack consistent theoretical 
justification. The higher-temperature peak (ii) shows many attributes of a thermally acti- 
vated relaxational process f^, but a calculation shows that to fit experiment, different sets of 
relaxational processes are needed for different z/ 0. Also the plateau region (iii) is difficult 
to explain by a thermally-activated relaxation process since numerical fits require unphys- 
ically large attempt frequencies [^]. Further, thermal relaxation processes give attenuation 
that increases more slowly than quadratic with increasing z/, contradicting (iv). Thermal 
vibrations have been overlooked as a sound-wave damping factor on grounds that vibrational 
modes would need unreasonably large Griineisen parameters (7 ~ 200 for vitreous silica 0) 
to account for the measured F. Until now, however, there has been no numerical study to 
test this argument. 

In this paper we examine the role thermal vibrations play in the sound attenuation in 
glasses. We will use the term "vibron" to refer to any quantized vibrational mode in a glass 
0. Our analysis is restricted to the region firin ^ 1 (the so called Akhiezer regime 0), 
where Tin is the inelastic lifetime or thermal equilibration time of a thermal vibron. We show 
that the unusually strong coupling (measured by Griineisen parameters 7) between sound 
waves and the low-frequency resonant modes explains the features (iii) through (v). As for 
the interpretation of (ii), our calculation shows that confusion arises because there actually 
are two different peaks. One is caused by relaxational processes (not addressed here) and 
dominates below 1 GHz and another is due to thermal vibrations and dominates at the 
lowest hypersonic frequencies. A double peak structure should be expected at intermediate 
frequencies. There is some indication for such structure in measurements on vitreous silica 
0. Our calculation is also a prediction: the existing measurements on amorphous Si |§] 
report F at too low frequencies (300 MHz) to see contributions of thermal vibrations. But 
even at higher frequencies (say, 30 GHz) one may expect traces of thermally activated peaks 
due to various defects. Recently discovered amorphous Si with 1 at. % H in which 
tunneling (and perhaps also relaxational) processes are inhibited would be excellent to test 
our results. 
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In the Akhiezer regime a sound wave passing through a sohd can be attenuated by 
two processes [|^. First, if the wave is longitudinal, periodic contractions and dilations in 
the solid induce a temperature wave via thermal expansion. Energy is dissipated by heat 
conduction between regions of different temperatures. Second, dissipation occurs as the 
gas of vibrons tries to reach an equilibrium characterized by a local (sound-wave induced) 
strain. This is the internal friction mechanism. To establish the relative importance of the 
two processes, consider order-of-magnitude formulas Th ~ {VL^ / pv^){KTa^ p^v"^ / C'^) and Fj ^ 
pf^) (CTr in7^) for the heat conductivity and internal friction processes, respectively |]10| . 



Here p is density, C specific heat per unit volume, v sound velocity, k thermal conductivity, 
and a is the coefficient of thermal expansion. The ratio TfjTi ^ (/tQ;^p^f^)/(C'^rin7^) 
becomes more intuitive when putting a ~ Cl/B {B ^ pv"^ is the bulk modulus) and 
K ~ CD where D is diffusivity. Then F/i/Fj D/{v^t\^). The factor f^rin measures the 
ability of vibrons to absorb energy from a sound wave of velocity v. The difference between 
a glass and a crystal lies in the values of D and Tin. In crystals D f» w^Tjn, that is, energy is 
carried by phonon wave packets with group velocity v. The ratio F/j/Fj is then of order unity. 
In glasses energy is transferred by diffusion (spreading rather than ballistic propagation of 
wave packets ||ll|) and D is not related to Tin One of the reasons the contribution to 



Fj of thermal vibrons was previously underestimated is that Tin was guessed from thermal 



conductivity p[; this gave too small t-^^. For amorphous Si D ~ 10 m /s [12], f 8 x 10 
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m/s, and Tin ~ lO^^^s give F/^/Fj 0.02. Since these are typical values, Vh can be 
neglected. This is consistent with experiment: compared with crystals glasses have smaller 
K and yet F can be larger p|. 
Internal 

friction leads to sound-wave energy attenuation [jlO| F = {^'^ / pv^q'^)riai3^sQa^i3(l'Y^s, where 



Va/B-yS is the internal friction tensor with cartesian coordinates a...S and q (e) is the wave 
vector (polarization) of the sound wave. Summation over repeated indexes is assumed. We 
will evaluate F for both longitudinal (L) and transverse (T) sound waves with wave vectors 
averaged over all directions: 

~ 1 K 3 iVaaf3f3 + '^VaPap), (l) 



15pvl 
30pw^ 



— Qn„„ 3 i'^VafiaP — VaaPfs)- (2) 



The coefficients rjap-ys are the real part of a complex tensor rja/s-ys which can be obtained by 
solving a kinetic equation in relaxation time approximation 

^ ^ 7a/3775 - (7^/37^5 + 7^/3775) /2 , . 

The summation is over all vibrational modes j; cj and tj denote mode specific heat and 
relaxation time. The Griineisen tensor —7^^ is the relative shift of mode frequency ujj per 
unit strain e^/j; 7 is the mode average of 7-' weighted with Cj/{1 — iQrj). The applicability of 



kinetic theory to the problem of internal friction was justified by DeVault and coworkers 
who obtained 77 from a microscopic theory as an autocorrelation function of the momentum 
current density operator. Remarkably, the microscopic theory shows that the momentum 
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current in a solid is not monopolized by ballistically propagating vibrational modes as in 
the case of the energy current. Nonpropagating (even localized) modes can contribute as 
much as propagating ones to the momentum current. One consequence is that the concept 



of "minimum kinetic coefficient," as introduced for electrical or heat [0 conductivity 
of disordered systems, is not realized for internal friction. We generalized DeVault's 
theory to include internal strain, the atomic rearrangements in a strained solid. We found 
that internal strain affects internal friction only by modifying 7-' , as in the case of thermal 
expansion |T^: 7-' now reflects the change between the initial mode frequency and the 



frequency of the mode after the rescaling (scaling parameter 1 + e) plus the rearranging of 
atomic positions (to achieve a new equilibrium at strain e). Internal strain is very important 
for thermal expansion of glasses we will show that it is important for r) (and F) as well. 

We calculate t] and F for the model of amorphous Si based on the Wooten- Winer- Weaire 
atomic coordinates and Stillinger- Weber interatomic forces |19|, with 1000 atoms ar- 



ranged in a cube of side 27.549 A with periodic boundary conditions. Diagonal Griineisen 
parameters 7^a/3 = {■yli + 722 + 733)/3 for this model were given in Ref. Jl^, transverse 
7q,/3 are calculated here. Vibrational lifetimes tj are extracted from their 216-atom version 
values (see also Ref. [^). The model has sound velocities vl = 7640 m/s and vt = 3670 
m/s|2g. 

Figure |l] shows the calculated F(z/) for longitudinal and transverse sound waves in amor- 
phous Si from 10 MHz to 1 THz at 300 K. The attenuation F ~ z/^ up to about 100 GHz, 
where the condition for the applicability of kinetic theory, Qr^^ ^ 1 reaches its limit (rin ^ 
1 ps). Our calculation is not valid beyond this point. In comparison, the measured attenu- 
ation of longitudinal waves in vitreous silica grows quadratically with 1/ up to at least 400 
GHz 0] suggesting that Tin in vitreous silica is several times smaller than in amorphous Si. 
This is not surprising since Si is remarkably harmonic: room temperature heat conductivity 
of crystalline Si is larger by an order of magnitude than that of quartz |2^, and a similar 
relation may hold for the corresponding Tin of the glassy phases. 

More surprising is the comparison with crystalline Si. Figure |l| shows that F^ is simi- 
lar for the amorphous and crystalline cases (measured F for vitreous silica is several times 
larger than for quartz 0). One would naively expect the sound attenuation in a glass to 
be much smaller than in the corresponding crystal since, owing to a distribution of bond 
lengths and bond angles, anharmonicity of the glass is higher (and Tin smaller). The same 
interatomic potential, for example, yields Tin for high-frequency phonons in crystalline Si at 
300 K about five times larger than in amorphous Si . The reason why F in glasses can be 
of the same order or even higher than in crystals is the internal-strain induced anomalously 
large Griineisen parameters of the resonant modes [|r^ (see also Fig.^. (Resonant modes are 
low-frequency extended modes whose amplitude is unusually large at a small, typically un- 
dercoordinated region |l3,|l.) Atomic rearrangements caused by internal strain are largest 



in the same regions of undercoordination where the resonant modes have largest amplitude 



17| . This leads to high sensitivity (measured by 7) of the frequencies of these modes to 
strain. If the internal strain is neglected, the sound attenuation is an order of magnitude 
smaller, as seen in Fig. |I|. (Since the resonant modes have low frequencies, their tj is longer 
than an average Tin; this adds even more weight to these modes.) Fewer than one percent of 
the modes are capable of increasing F by a decade! We believe the measured F for vitreous 
silica is also caused by the strong coupling of sound waves and resonant modes. Vitreous 
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silica is a much more open structure than amorphous Si so the number of resonant modes 
should be higher, bringing F above the crystalline value. 

Another interesting feature in Fig.|^ is the relative attenuation strength for longitudi- 
nal and transverse sound waves. While our model of amorphous Si gives Ti/Tt ~ 1/3 
at 300 K, the measured ratio for crystalline Si is reversed: Tl/Tt ~ 3 |^. This again 
shows how differently is sound attenuated in glasses and in crystals. The ratio Tl/Tt can 
be written as {vt IvlY{;-)\/'-^t)^ where 7l and 7t are effective Griineisen parameters. A 
crude way to estimate and 7!^, suggested by Eqs. |I] and ^ is to take mode averages of 
{iLlpp + 27a/37a^)/15 and (37^^7^^ - iLlpp)/^^- Our model gives il ^ 2 and 7^, ^ 1. 
The ratio T^/Tt is then about 1:3, in accord with the full calculation. Assuming the same 
ratio ~ 3 for vitreous silica (wl = 5800 m/s and vt = 3800 m/s), transverse and 

longitudinal waves are attenuated about equally. This is observed in experiment [Q]. The 
explanation of the measured Tl/Tx in crystalline Si can be found in Ref. . 

In Fig. 1^ we plot F(T) for different u. A remarkable feature is a peak at about 20 K 
at IMHz and below. As u increases, the peak shifts towards higher T and vanishes above 
4-5 GHz. Two factors cause the peak, (a) The sum J2jCj{'y^)'^ saturates at much lower 
temperatures (about 50 K) than the model Debye temperature To ^ 450 K ||2^. This is 
because the relevant j are resonant modes with small frequencies, (b) For low-frequency 
modes Ttj, after increasing linearly develops a peak, before going constant [much like T{T) 
itself]. As the temperature dependence of F follows J^jCji'y^)^ Ttj, the peak appears. At 
large u the peak vanishes because of the factor 1/(1 + i7^r^) in Eq. ^. At T above 100 K 
F(T) is nearly constant, as observed in experiment as a plateau (iii). This again follows 
from (a) and (b). 

We are not aware of any experiment with which we could compare our calculations. The 
measurement of F(T) of sputtered amorphous Si films reported in Ref. 0, for example, was 
performed at 300 MHz. This is too low to see any contributions from thermal vibrations. 
The whole temperature spectrum is dominated by a single peak of the type (ii), except at 
very low temperatures. This peak is expected to increase linearly with u, until thermal 
vibrations become relevant (roughly at 10 GHz), causing a plateau (iv) that increases as z/^ 
at higher frequencies. Even at smaller frequencies one may see some vibrational contribution 
to F(T) at large enough T, since the thermally activated peak decreases as 1/T at large T. 

Anomalous low T thermal expansion already suggested I^S] very large 7 values for low 
uj modes. Our large 7 values [0] agree nicely with trends in a{T). Like thermal expansion, 
F should be strongly sample and model dependent. There is evidence that our highly 
homogeneous model of amorphous Si becomes free of resonant modes when the number of 
atoms grows to infinity. That means an infinite model would predict F about a decade smaller 
than calculated here. Amorphous silicon, however, can be prepared only in thin films where 
voids and other inhomogeneities are unavoidable. Voids loosen the strict requirements of a 
tetrahedral random network (for example by introducing free boundary conditions). Then, as 
in our finite models, regions of undercoordinated atoms will allow the formation of resonant 
modes. While this issue for amorphous silicon will be ultimately settled by experiment, our 
calculation combined with the existing data on vitreous silica strongly suggests the reality 
of resonant modes. 

Our final note concerns the mode dependence of transverse Griineisen parameters like 
7i2. Similarly to volumetric 7aa/3 fT^, transverse 712 in Fig. |^ (713 and 723 look the same) 
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is unusually large for resonant modes and have scattered values for high-frequency localized 
modes. (More resonant modes have 712 negative than positive which suggests that resonant 
modes are trapped at highly anisotropic undercoordinated regions whose sizes change under 
shear [l^.) The 15 — 70 meV vibrons (diffusons 0) have 712 ~ (average magnitude 
0.02), while the corresponding 7ao/3 are of order unity ||T^. Such small values (zeros in an 
infinite model) are characteristic for diffusons, which are extended modes whose polarization 
directions (atomic displacements) point, in general, at random. There remains only a short- 
range correlation between polarization directions which determines the diffuson's frequency 
ujii. If a shear, say, ei2 is applied, Ud changes to a;^(ei2). Since long-range order in the 
diffuson polarization is absent, uj'^{ei2) ~ co'^(— ei2), and 712 which is a linear coefficient in 
the expansion of u'^ in eu must vanish. 

We thank J. L. Feldman for helpful discussions. The work was supported by NSF Grant 
No. DMR 9725037. J. F. acknowledges also support from the U.S. ONR. 
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FIG. 1. Log-log plot of the sound attenuation r(cm~'^) at 300 K, as a function of sound-wave 
frequency z^(Hz). Calculated T are represented by lines (IS-internal strain), experimental data by 
symbols. L (T) stand for longitudinal (transverse) sound waves. 



FIG. 2. Calculated sound attenuation r(r) for amorphous Si at different frequencies. The thin 
dashed lines are for longitudinal waves with the labeled u in GHz. Plotted are rescaled values F/u'^ 
for 1/ measured in GHz. 



FIG. 3. Calculated transverse Griineisen parameters 712 for amorphous Si as a function of 
vibron frequency. Above the vertical line (ftsTl meV) modes are localized. 
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